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Abstract. We study the mean number of encounters up to time t, E^ (t) , taking place 

in a subspace with dimension d* of a d-dimensional lattice, for N independent random 

walkers starting simultaneously from the same origin. En is first evaluated analytically 

f\\ in a continuum approximation and numerically through Monte Carlo simulations in 

"pi ■ one and two dimensions. Then we introduce the notion of the intersection of strongly 

r^ anisotropic fractals and use it to calculate the long-time behaviour of En- 
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1. Introduction 

Despite its long history, the study of random walks still remains an active field of 
research [1, 2, 3, 4]. For the discrete random walk, introduced by Polya [5], a quantity 
of interest for applications is the mean number S± (t) of distinct sites visited up to time 
t by a single random walker. It evidently grows as \fi when d — 1. The asymptotic 
behaviour for d > 1 was obtained by Dvoretzky and Erdos [6] with S\ (t) ~ t/ In t at the 
critical dimension d c — 2 and S\(t) ~ t for d > 2. Exact values of the amplitudes on 
different lattices in d = 3 and sub-dominant contributions were later calculated [7, 8]. 
Many other observables associated with the geometry of the support (i.e. the set of 
visited sites) of a random walk and their fluctuations have been studied more recently 
[9, 10]. 

The peculiarity of the dimension d = 2 was first noticed by Polya [5] who showed 
that, with probability 1, a random walker returns infinitely often to the origin when 
d < 1 and escape to infinity when d > 2. This property is linked to the value of the 
fractal dimension of the random walk, df = 2. 

Instead of considering a single random walk, one may generalise to the case of iV 
statistically independent random walks starting from the same point. Such studies were 
first concerned with the properties of first passage times [11, 12, 13, 14]. The study of 
the number of distinct sites visited by iV random walkers, S'jv(i), was initiated in [15] 
where asymptotic expressions for N large were obtained. In [16, 17] some corrections 
to [15] were given and sub-dominant contributions were evaluated. 

The mean value of another interesting observable, W^(t), the mean number of 
common sites visited by iV random walkers up to time t, was studied in [18] where the 
asymptotic behaviour at long time was determined. Quite recently, exact expressions for 
the probability distributions of the number of distinct sites and the number of common 
sites visited by N random walkers starting from the same origin have been obtained in 
d= 1 [19]. 

Following the work of Majumdar and Tamm [18], we have shown how the 
asymptotic behaviour of Wjv(i) can be simply obtained exploiting the notion of fractal 
intersection [20]. The rule governing the fractal dimension df(N) of the intersection 
of iV statistically independent and isotropic fractals with dimensions df, embedded in 
Euclidean space with dimension d, is the following [21] 



d f (N) = d- d f (N) , d f (N) = min [d, NdJ] , (1.1) 

where df = d — df is the co-dimension [22]. The dimension of the intersection vanishes 
at d — d c , the upper critical dimension, when the sum of the co-dimensions is equal 
to d [23]. The mean number of common sites visited, Wn (t), behaves as the mass 
M ~ R d t( N ) of the fractal resulting from the intersection of the N random walks. The 
typical size at time t is R ~ \/i and the dimension of the intersection follows from (1.1) 
with df = 2. 

In the present work we study the mean number of encounters up to time t, E^it), 
of N independent random walkers starting from a common origin at t — 0. The events 
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in space-time, corresponding to the meeting of the N random walkers, take place at 
the intersection of the N fractals associated with the support of the random walks [24] . 
Since, when considered in space-time, a random walk is a strongly anisotropic fractal 
object, the scaling behaviour of E^{t) at long time is governed by the fractal dimension 
of the intersection of N strongly anisotropic fractals. 

A fractal object is strongly anisotropic when its growth in different directions is 
governed by different exponents, R» ~ R z ± . The random walk in space-time, with a 
dynamical exponent z = 2, is one of the simplest examples. Such a behaviour is not 
uncommon, it is observed in different contexts among which one may mention directed 
systems [25], Lifshitz points [26, 27] and some quantum critical points [28], critical 
dynamics [29] and non-equilibrium phase transitions [30, 31, 32]. 

The paper is organised as follows. In section 1, E^{t) is calculated analytically 
in a continuum approximation for encounters taking place in a subspace of dimension 
d* of the rf-dimensional Euclidean space. In section 2, E^(T) is evaluated numerically 
through Monte Carlo simulations for discrete random walks in one and two dimensions. 
In section 3, the rule giving the dimension of the intersection of strongly anisotropic 
fractals is established, generalising (1.1). Then the notion of fractal intersection is used 
to evaluate the long-time behaviour of E^(t). In the final section, we formulate some 
remarks and present possible extensions. 

2. Mean number of encounters of iV random walkers 

We consider N random walkers performing discrete, statistically independent random 
walks on a rf-dimensional hyper-cubic lattice Cd with lattice parameter a. The N walks 
start simultaneously at t — from the origin of Cd- The position of a lattice site is 
given by f = Yli=i x i ^ where the Ui are mutually orthogonal unit vectors and Xi = nia 
(n,i e Z). The walks are also discrete in time, t — nr. 

We want to calculate the mean cumulated number of encounters of the N random 
walkers, occurring on a sub-lattice Cd* of Cd, containing the origin, with dimension 
d* [33] such that < d* < d. 

Let I^(f,t) be an indicator associated with a site at f for a given time t with the 
following values: /^(r, t) — 1 when the N random walkers visit the site at r*at the same 
time t and ijv(r, £) = otherwise. The number of encounters of the N independent 
random walkers up to time t, for a realisation of the N random walks, is given by 
Y^t'=o J2rec *lN(r,t'). The mean cumulated number of encounters of the N random 
walkers at time t, E^(t), is obtained by taking the average of this expression so that 

t t 

E N(t) = £ £ <W,f)> = E E p*W) , (2-i) 

t'=0 rdC d * t'=o fec d * 

where P^{f,t') is the probability to find the iV random walkers at r at time t'. For iV 
statistically independent random walks we have P/v(r, t') = [Pi(r,t r )] N . 
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The N random walkers starting simultaneously from the origin at t — 0, their first 
encounter takes place at t' — and the first sum can be split as follows: 

t 

Mt) = i+EE[ p '(^] ff ' ( 2 - 2 ) 

Since we are interested in the long-time behaviour of E^it) we use the continuum limit 
where both f and t are treated as continuous variables. In this approximation 

where .D is the diffusion coefficient. Thus when t^> t one obtains 

'df [ d d *rexp(-Nr 2 /4Dt') _ Nd 



(2.5) 



where to — r /2 will provide a cut-off when needed and 
Collecting these results leads to 

Nd-d* ft 

E M (t) « I H - / dt'(t') (d *- Nd)/2 (2 6) 

^ N[t) l + T ( AnD y N d-d*)/2 N d*/2 J to at[t > ■ V*) 

Let B^{t) denote the time integral in (2.6). It is given by 

f l+(d*~Nd)/2 _ ,l+{d*-Nd)/2 ,* „ 

B N (t) = In (pj , d = d e . (2.7) 

Thus the long-time behaviour of the integral changes at the critical dimension d c . It is 
governed by the constant cut-off when d > d c and by the time- dependent upper limit 
otherwise. The asymptotic behaviour of the mean number of encounters follows with 

E N (t) ~ t 1 **-"*" 2 , 0<d< d c (N, d*) = *±1 , 

~ ln(t/t ) , d = d c , 

~ const. , d > d c . (2.8) 

When d* = d the critical dimension, which is the solution of d c (N, d c ) = d c , is given by 

4(iV) = ^- I , d* = d. (2.9) 

Actually (2.6) remains valid for d = d* = where -Ejv(t) ~ t/r. Then the walkers 
cannot move and each time step gives an encounter. 

The behaviour of E^(t) in one and two dimensions, as obtained in (2.6-2.7), is 
shown in tables 1-2 for different values of N and d*. The time- independent coefficients 
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Table 1. Behaviour of the mean cumulated number of encounters Ej^{t) in d = 1. 
d= 1 N 1 2 3 4 

d* = ~ V* - ln(t/* ) "03 - Poa t_1/2 "04 - /?o4 * _1 
d* = 1 t + 1 ~y/i ~ ln(t/to) "14 - /?i4 £~ 1/2 

Table 2. Behaviour of the mean cumulated number of encounters E^it) in d = 2. 
d = 2 TV 1 2 3 



ef =0 


- ln(t/*o) 


alz-ffiit- 1 


"03 P03 * 


d* = 1 


~Vi 


«?2-^2*" 1/2 


"?3-A 2 3^ 3/2 


d* =2 


t + 1 


- ln(t/t ) 


"23 -^3* _1 



are written as oc^ N and /S^-tv- When N — 1, E^{t) gives the mean number of times 
a single walker visits the subspace of dimension d* up to time t. When d* = d 
the encounters are counted without restriction on their location. In particular, with 
d* = d and N — 1, the walker has an encounter with himself at each time step, hence 
E N (t) =t+l. 

3. Monte Carlo simulations 

The mean number of encounters of discrete random walkers has been studied numerically 
via Monte Carlo simulations for different values of N, d and d* in the different regimes 
around d c , for d = 1 and 2. 

The hyper-cubic lattice can be decomposed in two sub-lattices with J2 i=1 xi even 
on one sub-lattice and odd on the other. When the random steps are between nearest 
neighbour sites only, the random walkers starting from the origin at t — will stay 
on the even (odd) sub-lattice when t is even (odd). Due to this parity alternation 
the probability of encounter is artificially multiplied by a factor of 2 with respect to 
continuum results. 

In order to avoid this effect, we break the parity alternation by allowing the walkers 
to wait on the same lattice site with probability 1/2 or to jump with probability l/(4d) 
to one of the nearest neighbour sites. Thus we use the following single step probability 
density in d dimensions 

1 1 d 

Pif) = 9^) + TjX^^ -1-0 ^) + S(f-aUi)] . (3.1) 

i=i 

With a = t = 1 the diffusion coefficient is 

(r 2 ) 1 
D = ±-t = — 7 . (3.2) 

2dr Ad V ; 
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Figure 1. (Colour online) Comparison between the analytical results in the continuum 
approximation (lines) and the Monte Carlo data (symbols) for the mean cumulated 
number of encounters E^it) on Cd* of N random walkers starting from the origin of 
Cd at t — 0. The values of N and d* are such that d is below the critical dimension 
d c and En grows as \ft. When N = 1, En gives the mean number of visits of the 
sub-lattice Cd*- 



The mean number of encounters was evaluated by averaging over 10 6 samples so that 
the statistical errors are negligible. 

The Monte Carlo data for -Ejv(t) are compared to the analytical results of (2.6) 
and (2.7) for values of d* and N leading to the different regimes: power law growth 
below d c in figure 1, logarithmic growth at d c in figure 2 and power law approach to 
saturation above d c in figure 3. 

A systematic deviation between the analytical results (lines) and the Monte Carlo 
data (symbols) becomes apparent at, and mainly above d c where the difference is no 
longer negligible compared to the values of E N . It is clearly due to the failure of 
the continuum approximation at short time. In order to put it in evidence, the mean 
number of encounters at time t, e n (t), which is given by dE^/dt in the continuum 
approximation, is shown in figure 4 for values of d* and N leading to the saturation 
regime. The continuum results converge to the Monte Carlo data after only a few time 
steps. 



4. Intersection of strongly anisotropic fractals 

The random walk, when considered in space-time, is a strongly anisotropic fractal 
growing as t in the time direction and as \/i in the d transverse spatial directions. Thus 
counting the encounters of random walkers is a problem of intersection in space-time of 
strongly anisotropic fractals. 
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Let us consider a strongly anisotropic fractal in d + 1 dimensions. It has a fractal 
dimension dj and a typical size R in the d transverse directions while the corresponding 
values are dn and R\\ in the longitudinal (temporal) direction. The mass of the fractal 

is 



M ~ R a f 






i?ll = fl 2 



CLfW — — 
Z 



'711 



its volume is 



V ~ i?\Rii ~ i? 



d+z 



so that the average density scales as 

<p) M 



ydt—d—z 



V 



rs-/ £\^ J " r\J 



R~ d f 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



which defines the fractal co-dimension for a strongly anisotropic object 

d^=d + z-d f >0. 

Consider now iV statistically independent fractals with a common origin and the same 
values of df and z. The average density of their intersection, (pn), will scale as (p) so 
that 

-NdJ _ R -d7(N) 



(Pn) ~ (p) 



N 



Br 



(4.5) 



where df(N) is the fractal co-dimension of the intersection of the N fractals. Thus the 
fractal dimension of the intersection is given by 

d f (N) =d + z-d f (N). (4.6) 



i 1 1 1 1 1 1 1 r 
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Figure 2. (Colour online) As in figure 1 but the values of N and d* are such that 
d = d c . The growth of the mean cumulated number of encounters En is logarithmic 
in t. 
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Figure 3. (Colour online) As in figure 1 for values of iV and d* such that d > d c . 
The mean cumulated number of encounters En tends to a small constant value which 
is approached according to a power law, £ -1 / 2 for d = 1 and for d — 2 (squares), i _1 
for d = 2 (circles and diamonds). 
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Figure 4. (Colour online) Mean number of encounters ejv(i) on Cd* of N random 
walkers starting from the origin of Cd at t = 0. The results obtained in the continuum 
approximation (lines) quickly converge to the Monte Carlo data (symbols) . 
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Since the mass of the fractal intersection cannot decay with R, df(N) is bounded from 
above by d + z so that: 



d f (N) = min[d + z,Nd f } > 0. (4.7) 

The co-dimension of the fractal intersection is given by the sum of the co-dimensions of 
the fractals involved, when this sum is smaller than d + z. 

For N isotropic fractals with dimension df in D = d + 1 dimensions, z — 1, 
df = D — df and, according to (4.7), the co-dimension of the fractal intersection is 
df(N) = mm[D,Ndf] > 0, in agreement with (1.1) for the Euclidean dimension D. 

We are actually interested in that part of the fractal intersection which belongs 
to a subspace with dimension d* + 1 in space-time, i. e., its intersection with this 
subspace. According to the rule of co-dimension additivity we have to add the co- 
dimension d — d* = d + 1 — (d* + 1) to Ndf in ((4.7)). The total co-dimension is 
then 



d f (N,d*) = mm[d + z,Nd f + d-d*}>0, d f (N,d) = d f (N) (4.8) 

and the appropriate fractal dimension reads: 



d f (N,d*) = d + z-d f (N,d*). (4.9) 

The upper critical dimension d c (N, d*) is reached when the sum of the co-dimensions 



in df(N,d*) is equal to d + z, leading to df(N,d*) = 0. Then we have d c + z = 
N(d c + z — df) + d c — d* so that: 

d c (N,d*)=d f -z + ^-^. (4.10) 

When d* = d, according to (4.4) and (4.7), we have d c + z = N(d c + z — df) which gives: 
dc(N) = z + N(d,- Z ) (4n) 

Below the critical dimension d c (N, d*) the fractal growth of the intersection is 
governed by 

M*(R) ~ R d f( N > d *) ~ R z+d*-N(d+z-d f ) ( d < <q ( 412 ) 

according to (4.9). 

The logarithmic growth at the critical dimension can be obtained by integrating 
the local density of the fractal intersection Pn{t) ~ r~ df ( N ^ over the volume V*(R) ~ 
R d *R\\ ~ R d * +Z so that 

R r R 



a 



M*(R)~ PN (r)dV*(R)~ r -d f (N)+d*+z-i dr _ ^_ 13 ^ 



Alternatively, one may integrate the fractal density Pw(r) ~ r d f( N ' d *"> over the 
unrestricted volume V(R) ~ R d+Z which gives 



R r R 



a J a 



M*{R)~ p* N (r)dV(R)~ r ~d f (N,d*)+d+z-i dr _ ( 414 ^ 
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The two integrals lead to the same result when d < d c since df(N) — d* = df(N, d*) — d 
according to (4.7) and (4.8). At the critical dimension df(N, d*) = d + z and 

r R 

M*(R)~ r' 1 dr ~ ]n(R/a) (d = d c ) (4.15) 

J a 

Above the critical dimension df(N,d*) in (4.8) sticks to the value d + z so that 
d f (N,d*) = and 

M*(R) ~ const. (d > d c ) (4.16) 

asymptotically. 

For the random walk we have R(t) ~ y/i and R\\(t) = t. The mass of the fractal is 
the number of steps, M(t) = t. Thus, according to (4.1), df = z = 2. With these values 
of df and z the critical dimensions in (4.10) for d* < d and (4.11) for d* = d agree with 
the values given in (2.8) and (2.9). 

The mean cumulated number of encounters at time t of the N random walkers in 
Cd* is given by the mass of the fractal intersection M*(R) with R ~ \ft. Inserting 
this value of R(t) in (4.12), (4.15) and (4.16) one recovers the asymptotic behaviours of 
E N (t) in the different regimes (see (2.8)). 

5. Concluding remarks 

According to (2.7), the exponent governing the time dependence of En involves N, d 
and d* through the combination d* — Nd. Thus changing d* by Ad*, N by AN and d by 
Ad does not affect the asymptotic behaviour of En it) [34] when Ad* = NAd + dAN. 
This is the origin of the correspondences appearing in tables 1 and 2. For example, in d 
dimensions up to time t, the mean number of simultaneous visits of the origin (d* = 0) 
by N random walkers has the same time dependence as the mean number of encounters 
of iV + 1 random walkers anywhere (d* = d), a well-known fact for N = 1. 

Another remark is in order concerning the regime d > d c where En is asymptotically 
constant. This behaviour follows from (4.8) where the minimisation leads to df(N, d*) = 
d + z so that df(N, d*) = 0. If the minimisation is not taken into account, the fractal 
dimension in (4.9) takes on negative values when d > d c . Actually, this negative fractal 
dimension governs the approach to saturation in (2.7). 

The notion of fractal intersection introduced here for identical strongly anisotropic 
fractals can be extended to the case of fractals with different values of the anisotropy 
exponent z and/or the fractal dimension df. 

To conclude let us stress that the notion of fractal intersection used here could 
be efficient in other similar problems. Among these one may mention the asymptotic 
behaviour of Wn and En for Levy flights [35, 36]. 
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